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We present a phenomenological theory for phase turbulence (PT) in Rayleigh-Benard convection, 
based on the generalized Swift-Hohenberg model. We apply a Hartree-Fock approximation to PT 
and conjecture a scaling form for the structure factor S(k) with respect to the correlation length £2- 
We hence obtain analytical results for the time-averaged convective current J and the time-averaged 
■ mean square vorticity Q. We also define power-law behaviors such as J ~ e M , fl ~ e x and £2 ~ e _I/ , 

where e is the control parameter. We find from our theory that (x = 1, v > 1/2 and A = 2u + v. 
These predictions, together with the scaling conjecture for S(k), are confirmed by our numerical 
results. 
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I. INTRODUCTION 



CO Rayleigh-Benard convection (RBC) has long been a paradigm in the study of pattern formation This system 

consists of a thin horizontal layer of fluid heated from below. There are three dimensionless parameters to describe 
the system j3j. The Rayleigh number R = gad 3 AT '/nv is the control parameter, in which g is the gravitational 
acceleration, d the layer thickness, AT the temperature gradient across the layer, a the thermal expansion coefficient, 
k the thermal diffusivity and v the kinematic viscosity. Under the Boussinesq approximation, only the density of the 
fluid is temperature dependent. Then the Prandtl number a = v j n is all one needs to specify the fluid properties. 
The third parameter is the aspect ratio T = L/d where L is the horizontal size of the system. When the Rayleigh 
O , number R surpasses a critical value R c , the fluid bifurcates from a static conductive state to a convective state, in 
which the velocity profile u = (u±,u z ) and the temperature-deviation profile 9 form certain self-organized patterns. 
y—( I Those patterns also depend on the boundary conditions (b.c.) at the horizontal surfaces of the container. The two 
^ ■ most studied b.c. in the literature are the rigid-rigid b.c, under which the fluid cannot slip, and the free-free b.c, 
C*") ' under which the fluid does not experience any stress. 

The patterns and the corresponding stability domain in RBC have been studied extensively in the classical work of 
Busse and his coauthors QJ^| in the (R, er, k) space with k the wavenumber. For free-free boundaries, Zippelius and 
Siggia H and Busse and Bolton || found that the parallel roll state is unstable against the skewed- varicose instability 
immediately above onset if a < 0.543. Busse et al. J7j further investigated the dynamics involved and conjectured 
a direct transition from conduction to spatiotemporal chaos (STC). This spatiotemporally chaotic state is called 
"jS • phase turbulence (PT). Recently we reported a large scale (T = 60) numerical simulation of the three dimensional 
hydrodynamical equations for a — 0.5 under the free-free b.c. |^,^|. From that simulation, we confirmed the direct 
transition to PT above onset and studied various properties of it. The patterns we found have very complicated 
1 ' spatial and temporal dependences. 

In this paper we develop a theory for PT based on the generalized Swift-Hohenberg (GSH) model of RBC [|o| p^[ . 
This model is derived from the three-dimensional hydrodynamic equations, but is much simpler to study both nu- 
q I merically and analytically and is widely used in theoretical studies. Numerical solutions of this model or its modified 
versions have not only reproduced most patterns observed in experiments but also resembled experimental results 
relatively well [QJl^-|lJ] . But there are some shortcomings in the model : The stability boundary of the model 
does not coincide with that of hydrodynamics for both rigid-rigid and free-free boundaries; it induces an unphysical, 
short-ranged cross roll instability. Even so, owing to its simplicity and its qualitative resemblance to real systems, 
5^ , this model is very valuable in studying RBC. 

There are two coupled equations in two-dimcntional space r = (x,y) in the GSH model for two order parameters 
ip(r,t) and u) z (r, t), which are, respectively, related to the convective current and the vertical vorticity |17| . In this 
paper we present our analytical calculations of the time-averaged convective current J = A^ 1 J drip 2 (f,t) and the 
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time-averaged mean square vorticity ft = A 1 J drui 2 (r, t) in PT, where F(t) represents the time-average of F(t) and 
A is the area of the system. We carry out our calculations in Fourier space so the total number of modes ip(k,t) 
considered is infinite. We apply a Hartree-Fock approximation (HFA) to PT in which four-point correlation functions 
are approximated by products of two-point correlation functions. We also assume the time-averaged two-point correla- 
tion function C^(fci, £2) = ^(kijt^^jt) obeys a Gaussian distribution, i.e., C^(k\, ^2) = J {hi) exp[— \{k\ + fc2) 2 £ 2 >], 
where £^ is a length determining the correlation between different modes. This approximation is based on the ob- 
servation that the patterns in PT are disordered on long length scales and quasi-ordered on short length scales. 
This suggests that the correlation between different wavevectors k% and k^ is small unless these vectors are close 
to each other. This approximation seems physically intuitive for PT, in which there is a random superposition of 
rolls of different orientations. Given the above assumptions, we derive J and in terms of the time-averaged and 

azimuthally averaged structure factor S(k) = tp*(k,t)ip(k,t)/J. We further assume that S(k) obeys a scaling form 
kS(k) — ^2^F[{k — k max )^2\, in which £2 is the two-point correlation length, T{x) is the scaling function and k max 
is the peak position of kS(k). Applying this assumption, we obtain explicit formulas for both J and Q. More pre- 
cisely, we find that J = Joe — J^S,^ 2 ' which depends on unknown but experimentally measurable parameters, where 
e = (R — R c )/R c is the reduced control parameter and Jo and J^ are both known. We also find that £1 = ojJ 2 /^ 
where u> is related to the width of the scaling function T(x). Furthermore, by assuming power law behaviors such 
that J ~ e M , £2 ~ e ~ v and fl ~ e A , we predict from our theory that fi = 1, v > 1/2 and A = 2/i + v for PT. This 
prediction and the scaling assumption for S(k) have been verified by our numerical solutions. 

Our paper is organized as follows. In Sec. |ll[ we introduce the GSH model in Fourier space and derive the basic 
formulas governing the time-averaged convective current J and the time-averaged mean square vorticity fl for any 
pattern in RBC. In Sec. [II we introduce the HFA of our theory and use it to calculate explicitly J and fl for PT. The 



results are expressed in terms of the structure factor S(k). Sec. IV includes two parts. We first postulate the scaling 
form of the structure factor S(k) and expand both J and Q in the leading order of 1/^2- We then define power law 
behaviors for J, £2 and Q in PT and compare the results from our theory and our numerical work j^Q, which agree 
very well for both the exponents and the amplitudes. In the last section, we summarize our results and discuss some 
related issues. 



II. BASIC FORMULAS 



In the GSH model, the order parameter t^{r,t) satisfies ]10|-|l2|,|l 



to 



dt^ + U ■ Vip 



(£o/4fc c 2 )(V 2 



ip-N[ip]. 



(1) 



Here N[ip] is the nonlinear term to be specified soon and U (f) is the mean flow velocity given by U(r) = V£(f, t) x e z , 
in which 



dt - aV 



v 2 c 



V(V 2 ,//) x vv> 



(2) 



In the GSH equations, the reduced Rayleigh number e = (R/R c ) — 1 is the control parameter, in which R and R c are 
the Rayleigh number and its critical value at onset. The Prandtl number a parameterizes the fluid. While k c is the 
critical wavenumber at onset, other parameters model the properties of the system. More precisely, one takes [JL1 18 
R c = 27tt 4 /4, k c = tt/V2, r = 2(1 + cj- 1 )/^ 2 , $ = 8/3tt 2 , and g m = 6. 

It is easier to analyze the GSH equations theoretically in Fourier space than in real space. By convention, we define 
the Fourier transformation and its inverse transformation of an arbitrary function F(r) as F(k) = i J dre~ lk r F(r) 

and F(r) — 'Yl^F^e^ ' r , where A is the area of the system. Note that F*(k) = F(—k) for any real function F(r). 
It is easy to check that Eq. (0) can be rewritten in Fourier space as 

T d t ^(k) + V(k) = r(e; k)$(k) - Nik), (3) 

where V(r) = tqU ■ and 

r(e;k) = 6-e (k 2 -k 2 ) 2 /4kl (4) 
The nonlinear N(k) term has been evaluated at onset [O, 
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N(k) = ^ g(k- kW(k 2 Mk 3 )i>{k + k2 - k 3 ), (5) 

where the coupling constant g(cosa) is given in Ref. fTlfl with a the angle between k and k 2 . Rigorously speaking, 
the exact forms of Eqs. ([[]), (|) and (|) are derived near onset and deviations from them in real physical systems are 
possible for large enough e. But we disregard such complexity and take them as our model for further study. 

One may take an adiabatic approximation (dt = 0) in Eq. (|J) by neglecting the first term on the left-hand side. 
This term is small in comparison with the other terms, which can be verified by applying the same perturbation as 
that in phase dynamics JlS) ] . With this approximation, it now is easy to solve Eq. (|2|) for C(fc), which indicates that 
the mean-flow field is slaved by the ip(r,t) field. We are also interested in the vertical vorticity ui z (r) = — V 2 £(r). 
From Eq. (||) , it is straightforward to get that 

Co z {k) = k 2 C(k) = k 2 ), (6) 

where, with an exchange of index k 2 — > k — k 2 , 

f(k; k 2 ) = (fc 2 - 2fc • k 2 )(e z ■ k 2 x k). (7) 
Applying these results, one may easily evaluate the mean- flow contribution to Eq. (^), which is given by 

V(k) = v(k;k 2 ;k 3 )r(k 2 )^(k 3 )^(k + k 2 -h), (8) 

where 

iT T r ^ 9mTo [e z ■ k X (k 3 - k 2 )][e z ■ k 3 X k 2 ]{kl - k%) 
v(k;k 2 ;k 3 ) = — -= =— . (9) 

2(7 |fc 3 -fc 2 | 4 

Notice that the coupling constant v(k;k 2 ;k 3 ) is zero under two conditions: (1) If all k allowed in ip(k) point at 
one single direction, say k; or, (2) if all k lie on one single ring, say |fc| = k. For this reason, ordered states such 
as parallel rolls, hexagons, concentric rings, etc., do not have significant mean-flow couplings. Furthermore, the 
coupling constant v(k;k 2 ;k 3 ) seems to have a pole at k 2 = k 3 . Assume that k 3 = k 2 + q with q very small; then 
v ~ (e z • k x q)(e z ■ k 2 x q)(k 2 ■ q)/q A . A pole normally exists (since v ~ 1/q) unless q\\k, or q\\k 2 , or q _L k 2 . 
In this paper, we will mainly focus on two global quantities: One is the total convective current defined by 

J i t) = J [dr^(r,t) = 52j(k,t) with jQe,t)=$*(k,t)i>(k,t)\ (10) 

k 

the other is the mean square vorticity defined by 

fi(t) = j J dfcoUr, t)= + k 2 ;k 2 )f(ki + k 2 ; k 3 ) 

Notice that f{k\ + k 2 ; k 2 ) ~ (k\ ~ fe|)(e z • k 2 x k\). So Q(t) — if all the wavenumbers allowed in tp{k) point at one 
single direction ±fc or lie on one single ring \k\ = k. In other words, the vorticity must be generated by couplings 
between modes of different k and k. From Eq. (H), it is easy to derive that 



T Q d t J(k, t) = 2r(e; k)J{k, t) - ^ \g(k ■ k 2 ) + v(k; k 2 - k 3 ) 



k 2 .k-.i 



4>*(k)i>*{k 2 )iP(k 3 )^(k + k 2 ~ k 3 ) + c.c.J . (12) 

In principle, this is the equation determining the structure of the convective current J(fc, t) which, however, is beyond 
our present goal. Now applying the relations tp*{k) = i){—k) and v(k;k 2 ;k 3 ) = —v(k;k 3 ;k 2 ), and, exchanging the 
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summation indices fc — > —k, £2,3 - ► — ^2.3 for the g terms and k — > fc + &2 — fej ^2 &3 for the w terms, one obtains 
from the above equation and the definition of J(t) that 

± To d t J(t)= £)r(e;fc)J(M) 

- ^ g{h-hW{h)$%h)Ukz)${k + h-h). (13) 

fcl,fc2,fe3 

This is the equation determining the total convective current J(t). Notice that the v terms vanish from this equation, 
which can also be derived directly from Eq. (0) by converting the corresponding integral in Eq. ( jl(i| ) into a surface 
term. In general, the v couplings affect J(t) implicitly by modifying its structure J(fc,i) unless, of course, v = 0. 

For stationary states, the convective current and the mean square vorticity are time- independent. This, however, is 
no long true if the state is spatiotemporal chaotic. For a spatiotemporal chaotic state, these two quantities normally 
fluctuate in time around some well-defined averaged values: see Refs. HPI . While the fluctuations appear chaotic 
in time, they are relatively small in comparison with their averaged values. For simplicity, we only consider the two 
corresponding time-averaged quantities in our theory. We now introduce the time-average operator T defined by 
TF(t) ee Fit) = lim T ^ +00 i J T dtF(t). Applying T to Eq. (|l3|) yields 

J2r(e;k)J(k,t) 

k 

ki ,^2:^3^4 

In the next section, we show how, under a certain assumption, to calculate the time-averaged convective current of 
PT from this equation. The time-averaged mean square vorticity can be obtained with T acting on Eq. (]Tl"|). For 

simplicity, we denote from now on J(k) = J(k,t), J = J(t) and fl = O(i). 

Finally we introduce the time-averaged structure factor defined by S(k) = J(fc)/J with X^'S'(^) = 1, and the 
corresponding averages (F)j: = S(k)F(k). With this notation, the first term in Eq. (|l4| ) can be rewritten as 
(r(e))j:J. If the fc-dependence and the angular dependence in S(k) can be separated, then it is more convenient to 
define S{k) = (2 7 r) 2 A- 1 S'(fc)$(a) with J™ dkkS(k) = 1 and J™ da$(a) = 1, where a is the angle between fc and 
some reference direction. Here the discrete fc lattice has been converted into a continuous one. So a proper phase 
factor has been taken into account. Notice also that <&(7r + a) = $(a) since S(— k) = S(k). Now the corresponding 
averages with respect to S(k) and $(a) are defined as (F) k = f™ dk kS(k)F(k) and (F) a = da$(a)F(a). 

III. HARTREE-FOCK APPROXIMATION FOR PT 
A. Convective Current 

In order to calculate the time-averaged convective current and the time-averaged mean square vorticity from Eqs. 
( |l4| ) and (|TT|), it is obvious that more information about the corresponding state is needed. For this purpose, we notice 
that the instantaneous patterns in PT are irregular and random in space [||. So, if we write tp(k,t) = p(k, t)e^( fc '*) 
with both /5(fc,t) and 0(fc,i) real, the phase 4>(k,t) seems to have a rather complicated, irregular time-dependence. 
This suggests that two modes of different fc's are poorly time-correlated. Because of this, we think it is justified 
to adapt a Hartree-Fock Approximation (HFA) to PT in which a four-point correlation function is approximated by 
products of two-point correlation functions, i.e., 

cj, 4) (fci, fc 2 , fc 3 , fc 4 ) = ^(fci, t)i>(k 2 , t)i>(k 3 , t)j>(ki,t) 

~ CV(fci, fc2)CV(fc3, fc.4) + CV(fci, fc3)CV(fc 2 , ki) + C t i ) (ki,k 4: )C t f,(k2,k 3 ), (15) 

where we denote 



CV(fci, fc 2 ) = i/>(fci,i)?/>(fc 2 ,i). 



(16) 



Apparently C^,(k, —k) = J(k). 

Now it is essential to know the behavior of C^(ki, fe 2 ) in PT. We notice that the patterns in PT are disordered at 
large scales but quasi-ordered inside some smaller domains. From these, we speculate that the correlation between 
two modes of different fc's are poor if their fc's are quite different but good if their fc's are very close to each other. 
Hence, we postulate a Gaussian form for C^,{k\,k 2 ) in PT, i.e., 



C^kxM) = J(ki)exp[^-(k 1 +k 2 ) 2 ^}, (17) 



where £0 is a length determining the correlation between different modes. We expect £f to be large. 

We now calculate the convective current J from Eqs. (|l4|), (|l^) and (|l7|). One may rewrite Eq. ( |T4| ) in a continuous 
Fourier space as 

(r(e, fe))gJ = y dkxdk^dk^dki + k 2 -k 3 ~ fc 4 ) g(fci • fc 2 ) 

xCft-^-fe,!,^). (18) 
Inserting Eqs. ( |l5|) and ([I?]) into this equation, one finds after some algebra that 

J= ^MM^, (19) 

in addition to the conduction solution J = 0. Here S(k) is assumed to be azimuthally uniform, and gpr is defined as 

2 r 

g PT = gl— 1) -| — / dag(cosa). (20) 
7f Jo 

Using the explicit formula given in Ref. |Q] for free-free boundaries, one has that 

g PT = 0.855922 + 0.0458144CT- 1 + 0.0709326a-- 2 , (21) 

where a is the Prandtl number. 

If one takes the limit £^ — > +00, the two-point correlation function C0(fci,fc 2 ) in Eq. ( |l7| ) reduces to 
[A/2~K&)J{kx)5^ This, by applying the inverse Fourier transformation, leads to the translation invariance 

of the the two-point correlation function C(ri,f 2 ) = ^(riifyip^jt)/ J, i.e., C(ri,r 2 ) — C{f\ — f^). Under the same 
limit, one also has £^ = [A/27:] 1 / 2 since C^,(k, —k) = J(fc). In this case, one has 



J= mA. (22) 

gpr 



From now on, we simply choose £^ = p4/27r] 



1/2 



B. Mean Square Vorticity 

We now use the HFA to calculate the time-averaged mean square vorticity. From Eq. (|TT|), one finds that 

A 1 f - - - - (2-k) 2 ^^^^ ^ ^ ^ ^ ^ ^ 
ft = T^js y dk 1 dk 2 dk 3 dk 4 , ^ 5{ki + fc 2 - fc 3 - fe 4 )/(fei + &25 k 2 )f(ki + k 2 ; k 3 ) 

xC^i{~k 1 ,-k 2 .k 3 ,k i ) 1 (23) 

where, from Eq. (0), 



/(fci + fc 2 ; k 2 )f(kx + fc 2 ; k 3 )8(ki + k 2 - k 3 - fc 4 ) 
^ ffm ~ fc l)( fe I - kj)(e z ■ ki x fc 2 )(e z • fc 3 x 4) 

~ 4(72 |fcl + k| a |fc 3 + fc*| 2 



<J(fci + fc 2 -fc 3 -fc 4 ). (24) 
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Notice that there is a singular point k± + k% = in the above expression: its value depends on how the point is 
approached. 

Following the same calculation as for J and assuming that S(k) is azimuthally uniform, one finds after some algebra 
that 

= ^- J dhhSih) ydfc 2 fc 2 5(fc 2 )A(fc i; fc 2 ), (25) 
where we have chosen £^ = [A/271"] 1 / 2 , and 

A(fc i; fa) = i|fc 2 - fe|| [fc 2 + fe| - |fc 2 - fe||], (26) 
which has a second-order singularity at fci = fe 2 and is due to the singularity in Eq. (p4|). 

IV. SCALING RELATIONS IN PT 
A. General 

To evaluate the convective current J and the mean square vorticity, f2, one must know the structure factor S(k) 
which, however, is beyond our present theory. We thus turn to phenomenological arguments. We define a two-point 
correlation length as 

t2=[(k 2 ) k -(k)l}- 1/2 . (27) 
Then we assume that the structure factor satisfies the following scaling form 

kS(k) =&F[(fc-km«,)6], (28) 

where k max is the peak position of kS(k) and T(x) is the scaling function satisfying f_ dx !F{x) — 1. [Since fc > 
in fcS'(fc), the lower limit for !F(x) is —k max ^2, which we approximate by — oo.] Inserting fc = k max + x^ 1 and Eq. 



(pq) into Eq. (27), one gets that (x 2 ) x — (x) 2 = 1, where we have used the notation (F(x)) x = J_ QC dx T(x)F(x). It 
is also easy to see that (fc)fc = k max + £j (a;) x . 

For very large £ 2 , one may take fc = k max + x^ 1 in Eqs. (p2|) and (|25|) and expand the results in order of l/£ 2 . It 



— — — ^ a j . S a Z ^ — M y \ j 

is easy to find from Eqs. (|j) and (|22|), and from k m ax = k c + O(e) that 

2 



J = 

9pt 



2 Co 
(X )x-p2 

S2 



+ 0(l/^ 3 ,e/6,e 2 ), (29) 



for small enough e and large enough £ 2 . This expression depends on two unknown but experimentally measurable 
quantities £2 and {x 2 ) x = 1 + (.t) 2 . If T{x) is symmetric, then (x) x — and (a; 2 )^ = 1. In general, however, one has 
(x 2 ) x > 1. Notice also that although mean-flow couplings are not explicitly present in Eq. (|^), they affect the value 
of the convective current via the structure factor S(k) [see Eq. (|l2|)] and the two-point correlation length £ 2 . For PT, 
since £ 2 ~ §£o/ e5 (H, the contribution from £ 2 reduces the value of J quite significantly. This feature is absent in the 
simple theory presented earlier || . 

One may evaluate the mean square vorticity in the same way. From Eqs. (p5|) and (|26|), one gets for free- free 
boundaries that 

n^ 9 ^{\ Xl - X2 \) x ^^, (30) 

where we have used k max w fc c . The quantity (|xi — X2\) Xl>X2 is related to the width of the scaling function T{x) 
and, since (\x\ — x 2 |) Xl ^ 2 < yj JJxl — X2) 2 )x u x 2 = we expect (\xi — X2\) Xl ,x 2 to be of order of unity. Clearly, by 
phenomenological arguments, we can express J and O in PT in terms of measurable quantities. 
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B. Power Laws 



For PT, we further assume power law behaviors for the two-point correlation length, the convective current and 
the mean square vorticity such as 



6«6,oe _!/ , J~Jo^ and fi » f> e A - 
Then, from Eq. (|3(]), we find the following scaling relation 

A = 2[i + v. 

Recalling Eq. (p3), one obtains that 



2 

gpT 



P 2 

SO ,2;/ 
S2,0 



Since J is positive by definition, the values of the exponents satisfy 

fi = l, v>\/2 and A = 2 + i/>5/2. 



It is very likely that v = 1/2, hence, A = 5/2. If so, then one finds from Eqs. (29) and (pOD that 



Jo = 

9pt 



2.0 



„2 1.3 j2 



and CI 



o 



(|xi - x 2 |) 3 



(31) 



(32) 



(33) 



(34) 



(35) 



which depend on three phcnomcnological parameters £2,0^ {x 2 ) x and (\x\ — X2\} Xl ,x 2 - If f > 1/2, then Jo = %/9pt 
since the t 2v term in Eq. ( |33| ) contributes only to the leading correction to scaling. It is interesting to notice that the 
amplitude equations coupled with mean-flow jL2j predicts for free-free boundaries that fl ~ e 5 / 2 for almost perfect 
parallel rolls. 

We now verify our predictions for the power laws in PT by our numerical solutions. We have carried out large-scale 
numerical calculations of the three-dimensional Boussinesq equations under free- free boundaries for fluids of a = 0.5 

f,|). We have confirmed in Ref. M that the structure factor in PT satisfies the scaling form (|2^). From Table 
one can see that our theoretical and our numerical results are in very good agreement for the exponents. The 
scaling relation Eq. (B2J) is confirmed within our numerical uncertainties. The comparison between the corresponding 



amplitudes, however, is only moderately successful. Calculations of £2,01 ( x )x an d £2!) a 



/xi,£C2 are obviously beyond 

the present theory, so we take our numerical result for £ 2 ,o- Since our numerical results for (x 2 ) x and (\xi — X2\) xi .x 2 
are too sensitive to the large value cutoff to be meaningful, we assume equalities in (x 2 ) a - 1 1 

35|) and (pl|), one gets Jo 



(\xi - x 2 \) 3 



< 



y/dxi - £ 2 ) 2 ) 3 



V2. From Eqs. 



1 + (x) 2 x > 1 and 
0.972, which is about 20% 



larger than the numerical value. A non-zero value of {x) x will apparently reduce the theoretical value of Jo in the 
right direction. On the other hand, one finds that Slo — 454.7(|aci — £2|)xi,x 2 — 643.0 as an upper bound, which 
is about ten times larger than our numerical result. Nevertheless, we note that while our theory is based on the 
two-dimensional GSH equations, our numerical calculations are done for the three-dimensional Boussinesq equations. 
Although the former is very good in reproducing qualitative features of RBC, it may not be quantitatively accurate 
111 modeling RBC §. So one should be cautious in comparing the results from the GSH equations with those from 
real experiments or those from numerical calculations with hydrodynamical equations. 



V. DISCUSSION 



Our phenomenological theory for PT in RBC depends on two basic assumptions. In Sec. Ill, we adapt a Hartree- 
Fock approximation to PT and postulate a Gaussian form for the time-averaged two-point correlation function. In 
Sec. IV, we further assume that the structure factor satisfies a scaling form such as kS(k) = ^^[{k — fc ma i)6]- I n 
comparison with similar scaling forms in critical phenomena, critical dynamics and phase ordering, we find it necessary 
to replace k with k — k max in the scaling form. The physical origin of this replacement is due to the fact that patterns 
in RBC have an intrinsic wavenumber, which is close to k c . By the same reason, we find it necessary to seek the 
scaling form of kS(k) instead of S(k), where the k factor comes from dk = kdkda in two-dimensional fc-space. The 
physical origin of this scaling invariance in PT is yet unknown. In Sec. IV, we have confirmed the scaling form of 
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S(k) within our numerical accuracy. It is not clear currently whether the scaling form breaks down for small or large 
k. 

In summary, we present a phenomenological theory for PT in RBC. We calculate analytically the time-averaged 
convective current J and the time-averaged mean square vorticity O as functions of e and £2- We believe that our 
theoretical results will be useful in understanding the complicated behavior of STC in RBC. We also believe that 
our theory provides a new approach to STC and also raises some interesting questions. For example, how can one 
calculate the structure factor S(k) and the two-point correlation length £2 analytically? Is it possible that certain global 
quantities in STC form a complete set in the same way as temperature, pressure and density do for thermodynamic 
systems? Can we derive some effective variational principle in terms of global quantities? How far can we apply the 
ideas in critical phenomena to study STC? Since our assumptions are quite general, it will also be interesting to see 
whether our theory can be generalized to STC in other systems Jj|]. 
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A* 


V 


A 


£2,0 


Jo 


n 


Numerics 


1.034 ±0.025 


0.472 ± 0.016 


2.55 ±0.10 


0.82 ± 0.04 


0.787 ±0.019 


70.1 ± 1.0 


Theory 


1 


1/2 


5/2 




0.972 


643.0 
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